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Review of Cosmic Ray experiments with underground detectors 

E. Scapparone'^ 

^INFN, LNGS S.S.17 km 18+910, 61070, Assergi (AQ), Italy 

The most important underground detectors addressing Cosmic Ray physics are described, with a special em- 
phasis on the description of the used technology. 



1. Introduction 

The first generation of underground experi- 
ments ran during the '80s, triggered by the GUT 
theories. The most important results they ob- 
tained were the 1MB analysis on proton de- 
cay, ruling out the minimal SU(5) model, the 
^B solar neutrino observation by Kamiokande 
and the Supernova SN1987 detection. The im- 
pact on the technology development was rele- 
vant. Let me quote for instance the NUSEX 
experiment [y : this detector proved the relia- 
bility of digital calorimeters based on plastic 
streamer tubes, opening the road to the construc- 
tion of the hadronic calorimeters for the LEP 
experiments(ALEPH|l, 0PAL|1). 1MB and 
Kamiokande showed that the water Cherenkov 
technique was robust while LSD pointed out the 
possibility of using liquid scintillator to search 
for Supernova explosion. The main limitation 
of these detectors was the collected statistics . 
The natural development was therefore the con- 
struction of similar detectors enlarged by an or- 
der of magnitude in volume or in area. The ex- 
perience gained with NUSEX allowed the con- 
struction of MACRO [Q, with an acceptance in- 
creased to more than two orders of magnitude. 
The Kamiokande and the 1MB Collaborations 
joined Super-Kamiokande, reaching a mass of 50 
Ktons and LSD inspired LVD at Gran Sasso, ac- 
tually running with 600 tons of liquid scintillator. 
The most important questions are: how much 
the technology of these experiments is advanced 
? Do the most important technological limita- 
tions come from the large volume/ area used or 
this kind of physics requires just a modest detec- 
tor performance ? Could they benefit of more 



sophisticated techniques ? I will try to answer to 
these questions reviewing the main underground 
detectors for Cosmic Ray physics actually run- 
ning. 

2. MACRO 

Let me start with the MACRO experiment, lo- 
cated in the HALL B of Laboratori Nazionali del 
Gran Sasso. This is a multi purpose experiment 
for GUT magnetic monopole search, atmospheric 
and Supernova neutrino and Cosmic Ray stud- 
ies. The detector consists of 14 horizontal lay- 
ers of plastic streamer tubes(PST) and 3 hori- 
zontal layers of liquid scintillator (LS). Each side 
is closed with a sandwich of 6 PST and a layer 
of LS. The cell size of the streamer tube system 
has a cross section of 3x3 ctti^, larger than the 
NUSEX one (1x1 cm^). The readout was per- 
formed in NUSEX by using orthogonal pick-up 
strips, while in MACRO the wires and a layer 
of pick-up strips, forming an angle of 26.5° with 
the direction orthogonal to the wires, are read. 
The larger cell size, due to the dominant effect of 
the multiple scattering of the muons through the 
rock, doesn't affect significantly the pointing ac- 
curacy. Fig.l shows the angle between multiple 
muons detected by MACRO. Since high energy 
multiple muons can be considered parallel within 
few milliradians, the distribution is the folding of 
the muon multiple scattering with the detector 
angular resolution. The 68% of the events has 
an angle 9<l.r and therefore cr0=l.l7%/2~O.8°. 
The distribution is largely dominated by the mul- 
tiple scattering, whose contribution at this depth 
is around CTg^^~ 0.6°. This is a typical exam- 
ple in which the granularity of the detector is 




Figure 1. Space angle A9 (degrees) between the di- 
rection of tracks in two muon events. The dashed 
area contains the 68% of the distribution. 

limited by an external source of indetermination; 
an improvement of the detector space resolution 
would result just in a small improvement of the 
muon position measurement and pointing accu- 
racy. Fig. 2Q shows the MACRO magnetic 
monopole flux upper limit. Different techniques 
are used to identify a monopole in the various 
/3 regions explored. At (10-'*</3<10-3) the liq- 
uid scintillator informations are very helpful: the 
magnetic monopole signature is a train of pho- 
toelectrons, with a duration proportional to the 
monopole velocity. This is a delicate search, since 
background contamination could mimic the sig- 
nal: a simple ADC/TDC system is not enough! 
MACRO takes advantage of a Waveform Digitizer 
System, with a sampling frequency of 200 MHz. 
The system is calibrated using LED light. Fig. 3 
(l2| shows a muon crossing three layers of liquid 
scintillator and the respective Waveforms. The 
study of upgoing muons, originating from neu- 
trino interaction below the apparatus, requires 
high detector performance too. The discrimina- 
tion against the much more abundant downgo- 
ing muons is obtained by using the time of flight 
system (TOF). The rejection factor required for 
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Figure 2. MACRO monopole flux upper limit. 

this analysis is ~10~"'. Fig. 4 shows the 1//3 
distribution obtained by MACRO. The left peak 
conies from upgoing muons, while the right peak 
conies from downgoing muons. The TOF reso- 
lution is better than Ins, allowing a clear sepa- 
ration between the two peaks. Of course higher 
timing performance with TOF systems have been 
reached, for instance a resolution crt~100 ps has 
been obtained by the AMS TOF system ||l||. Nev- 
ertheless, it must be stressed that in this case 
a (Tt<lns resolution is obtained on a huge scin- 
tillator mass (600 tons), making use of 12 me- 
ters length counters, running for a long period, 
more than 8 years. As far as the physics results 
are concerned, such TOF system performance al- 
lowed the study of atmospheric neutrino anomaly. 
MACRO pointed out iQ a discrepancy between 
the real data and the Monte Carlo prediction 
for atmospheric i/^ flux . The ratio between the 
data and the Monte Carlo prediction for upgoing 
muon flux was R=0.74±0.036(stat)±0.046(sys)± 
0.13(theoretical). An improvement of the agree- 
ment between Monte Carlo and the MACRO 
data is found supposing an oscillation Vfi^>Vr 
or v^-^Vs with Ar7i^=2.5-10~^eF^ and sin^ = 1.0. 
Further informations about neutrino oscillation 
could come from the measurement of the muon 
energy. Such a measurement for through-going 
muons is very difficult with the present genera- 



Figure 3. Muon signal in the MACRO scintillator 
obtained with the 200 MHz Waveform digitizer . 
tion of underground experiment. This is a typi- 
cal situation in which a more advanced technol- 
ogy would have been useful. The average energy 
of upgoing muon is ~ 20 GeV, while the average 
energy of downgoing muons at LNGS rock depth 
is Ef^~300 GeV. A first attempt to measure the 
muon energy is due to the NUSEX and to the 
LVD Collaborations. 

3. LVD 

LVD is a multipurpose detector consisting of 
a large volume of liquid scintillator interleaved 
with limited streamer tubes in a compact geom- 
etry [^. The apparatus has a modulare struc- 
ture that consists of aligned towers of 38 modules 
each. Every module contains 8 liquid scintilla- 
tor counters of dimension 1.5x1x1 m^ seen by 
three photomultipliers. Taking advantage of the 
L shaped detectors used as tracking system, LVD 
is able to measure the depth-intensity curve up 
to 20 Kmwe( Kilometers of water equivalent). As 
far as the measurement of muon energy is con- 
cerned, LVD measured the quantity <AE/AL> 
as a function of the rock depth h , where L is the 
track length in the scintillator counter and AE is 
the energy lost in the counters. For h<8 Kmwe, 
this observable is expected to increase with h, 




Figure 4. TOF measurement obtained by 

MACRO: right peak downgoing muons, left peak 
upgoing muons. 

since a selection of higher rock depth corresponds 
to higher threshold for Cosmic Ray muon energy. 
For h>8 Kmwe, the muon flux coming from at- 
mospheric muon becomes negligible. Observed 
muons come from atmospheric neutrino interac- 
tion in the rock in the neighborhood of the de- 
tector: Vfj^+N—^fi+X. The different origin of the 
muons, manifests for instance in the change of 
their average energy. This is small signal, dif- 
ficult to observe. Fig. 5 shows <AE/AL> as a 
function of h. Although the statical error for high 
depth is large, at h>8 Kmwe a clear decrease of 
<AE/AL> is visible in the LVD data. 

4. Precision measurements 

Cosmic ray detectors are usually supposed to 
perform measurements only with 10-20% accu- 
racy. I would like to show two examples of 
small signals observed deep underground. Let 
me start with the measurement of the distri- 
bution of the distance between muon pairs de- 
tected deep underground, the so called "deco- 
herence" . This distribution depends on differ- 
ent cosmic rays features: the C.R. cross section, 
the muon parent mesons pt distribution, the mul- 
tiple scattering of the muons through the rock. 
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Figure 5. Measurement of <AE/AL> in the LVD 
scintillator as a function of the slant depth. 

A special care has been devoted to the study 
of the tail at high distance of this distribution, 
to search for anomalous high Pf. On the con- 
trary the surprise came from the low distance 
region. Fig. 6 shows the decoherence function 
observed by MACRO. A good agreement is evi- 
dent between the MACRO data and the HEMAS 
Monte Carlo expectation. Nevertheless, a dis- 
crepancy of 34% is found for muon separation 
D<80 cm. The same effect has been pointed 
out by the LVD Collaborationllq]: again the real 
data shows an excess respect to the Monte Carlo. 
An explanation of this effect has been proposed 
by the LVD Collaboration: muons at short dis- 
tance can come from the direct muon pairs pro- 
duction by muons: fi+'N^fi+'N+fj.^+fj,~ . It has 
been pointed out|^ that such process is sup- 
pressed, respect to the electron pair production, 
by a factor ml/rni only for small \=E^/ E^. At 
high V the cross section is not suppressed by 
such factor. Since the produced muon pair lies 
at low distance from the original muon, the re- 
sults is an excess of events at short distance in 
the decoherence function. A data analysis from 
MACRO fl^ shows that after the inclusion of this 
process in the muon transport simulation code, 
the agreement between data and Monte Carlo 
is restored. Another example of small signal 
pointed out in cosmic rays physics is the detection 
of photonuclcar interaction of muons deep un- 
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Figure 6. The low distance region of the MACRO 
experimental decoherence before and after the sub- 
traction of the i-i+N-^ii+N+ii'^ +fi^ process. 

derground. MACRO looked for charged hadrons 
with E>300 MeV, produced in the rock surround- 
ing the detector ||l^. A data analysis based on 
11,000 hours of live time, they found 1,938 candi- 
date events over a total sample of 9,544,318 muon 
events. The result was expressed in terms of the 
ratio Ri_i+h of the selected /x-1-hadrons events to 
the number of muon events in the same live time. 
After the subtraction of the background, they 
found for R^+u in real data and in the MC sim- 
ulations : 
- R°th^ = (1-91 ± 0-05sta ± 0.03,^,) • 10-4, 



R 



FLUKA 



(1.89±0.16,ta±0.02,j,,)-10- 



- R';fh''^ = (1-31 ± 0-14.ta ± 0.02,,,) • 10-5, 
confirming thus the obsolete treatment of the 
muon photonuclear interaction in GEANT3.15- 
3.21 pg|l pointed out inpOf. 



5. Super-Kamiokande 

The Super-Kamiokande experiment is the un- 
derground detector with the largest mass actually 
running. Events produced in 50 Ktons of water 
are reconstructed by 11,146 PMTs. The main 
goals of the experiment are the study of ^B so- 
lar neutrino, the search for proton decay and the 
study of atmospheric neutrino. As far as the at- 



mospheric neutrino is concerned, several measure- 
ments can be performed by SuperK: the u^/ve ra- 
tio, the measurement of stopping/through-going 
muons and the measurement of the v^ flux as a 
function of L/E, where L is the distance trav- 
elled by neutrino and E is the neutrino en- 
ergy. The analysis of the ratio R=(f,i/t'e)data/ 
{v^/ve)MC has been performed by several experi- 
ments. The NUSEX measurement [^ was affected 
by a large statistical error: R=(0.96t'.28)' while 
the FREJUS Collaboration reported a substan- 
tial agreement between data and Monte Carlo [nsl: 
R=(1.00±0.15(stat)±0.08(sys)). On the contrary 
a deficit of atmospheric v^ has been measured by 
two Cherenkov detectors. IMBm found a ratio 
R=(0.54±0.05±0.11) and KamiokandeQ mea- 
sured R=(0.60±0.05±0.05). The apparent mis- 
match between calorimetric and Cherenkov re- 
sults, was recently restored by the Soudan II 
calorimeter. The Soudan II Collaboration re- 
ported [§ a ratio R=(0.64±0.11(stat)j:°°^(sys)) 
based on an exposure of 3.89 kton-year, the 
largest for a calorimetric experiment. The 
analysis of SuperK was based on 33 Kton- 
year exposure p3|. They selected 4353 fully 
contained events and 301 partially contained. 
Only single ring events were used up to 
now. They splitted the events in two cate- 
gories: sub Gev events ( £'t,is<l-33 GeV) and 
multi-GeV events (£;^,s>l-33 GeV). They found 
R=(0.63±0.03(stat)±0.05(sys)) for the sub-GcV 
sample and R=(0.65±0.05(stat)±0.08(sys)) for 
the multi-GeV sample. These results, based on 
a very large statistical sample, confirmed the 
previous measurements. Nevertheless, the inter- 
pretation in terms of neutrino oscillation yields 
for Am^ a value 5•10-''eT/2<Am2<6•10-3eV"^ 
while the best fit for the previous experiments 
gave a larger Atti^, for instance Soudan II found 
Am^=l.l-10~^eT^^. The most promising way to 
study neutrino oscillation using atmospheric neu- 
trino seems to be the evaluation of i^ flux as 
a function of L/E: $^^= ^^.^ (L/E) ||. Fig.| 
shows the ($z.^)data/(*i/j,)MC as obtained by Su- 
perK. The agreement between the real data (cir- 
cles) and the Monte Carlo expectation(dotted 
line), obtained supposing a v^^Vt or v^^fVs os- 
cillation with parameters Am^=2.2-10~'^eF^ and 
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Figure 7. (^„^)data/('^ujMC as a function of 
L/E obtained by SuperK. 

s\n?{29)—l., is impressive. From a technical point 
of view, the SuperK analysis is based on single 
ring events only, therefore the reconstructed en- 
ergy E of the event relies only on the muon energy. 
The precision in the measurement of L is also re- 
duced, since the event kinematic is not completely 
defined. The most important goal of the next gen- 
eration of atmospheric neutrino detector is the 
improvement of the L/E ratio measurement. Dif- 
ferent ideas have been proposed. The measure- 
ment of the hadrons produced in neutrino events 
with a highly segmented calorimeter with mass 
~10Ktons or the use of a massive calorimeter 
(M> 30 kton) to measure fully contained events 
with any track multiplicity. 

6. Conclusion 

Usually the underground detectors for Cosmic 
Ray studies are much more coarse than acceler- 
ator experiments. Such a low granularity comes 
from the high area/volume required, imposing an 
upper limit on the detector cost and on the num- 
ber of electronic channels to be taken under con- 
trol. Moreover sometimes the coarseness of the 
experiment is fixed by the working conditions. 
Nevertheless, a good technical performance is re- 
quired to these experiments: 



- Long data taking. Underground Observatory 
looking for Supernova neutrino (LVD, MACRO, 
SuperK) have to run for a fong time (>10 years), 
avoiding any aging effect. 

- Flexibility. An high flexibility is required to 
these detectors: the physics to be investigated 
spans in energy over more than 3 orders of mag- 
nitude. At E~10MeV they look for Supernova 
neutrino while low energy hadrons produced in 
photonuclear interaction of muon have an energy 
of few hundreds of MeV. These detectors mea- 
sure atmospheric neutrino with E~l GeV, while 
the average energy of cosmic ray muons detected 
at Gran Sasso depth is £';i~ 300 GeV. 

- Precision measurements. In some analysis 
an high rejection factor against background is re- 
quired and signals as small as few 10^** are suc- 
cessfully measured. 

About the future the situation is quite different 
with respect to the '80s, when the most impor- 
tant lack of underground experiments was the 
poor statistics collected. The step forward was 
the construction of larger detectors, while now 
it is not possible to enlarge the underground de- 
tectors area/volume by an order of magnitude. 
Such mass increase is feasible only with under- 
water Cherenkov detectors. Moreover the under- 
ground space is limited and the safety rules be- 
came much more restrictive. Special care has to 
be devoted in the choice of the used materials. 
As far as the Cosmic Ray composition study with 
underground experiments is concerned, the inter- 
pretation of the results obtained by many experi- 
ments is nowadays much more urgent than the re- 
quest for further experimental data. The lack of 
accelerator measurements at high energy and at 
high rapidity region, introduces large uncertain- 
ties in the Monte Carlo models for the Cosmic 
Ray shower development. The trend for the next 
generation of underground detectors seems there- 
fore to be the construction of atmospheric neu- 
trino detectors with reasonable mass and with an 
increased granularity. The ICARUS detector [ll6| 
for instance is foreseen to run within the 1999 
at Gran Sasso with 600 tons. It will take ad- 
vantage of the superb space resolution to per- 
form an unambiguous analysis of the i-'^/i^e ratio. 
Other detectors have been proposed to improve 



the L/E measurements for atmospheric neutrino. 
Together with the Long Base Line experiments, 
they will hopefully give an answer on the nature 
of the neutrino, using more advanced technology. 
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